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Geometric properties of the set of quantum entangled states are investigated. We propose an 
explicit method to compute the dimension of local orbits for any mixed state of the general K x M 
problem and characterize the set of effectively different states (which cannot be related by local 
transformations). Thus we generalize earlier results obtained for the simplest 2x2 system, which 
lead to a stratification of the 6D set of TV = 4 pure states. We define the concept of absolutely 
. separable states, for which all globally equivalent states are separable. 

(N 
O 

Q 

I. INTRODUCTION 

> 

Recent developments in quantum cryptography and quantum computing evoke interest in the properties of quantum 
entanglement. Due to recent works by Peres Q] and Horodeccy there exist a simple criterion allowing one to judge, 
whether a given density matrix p, representing a2x2or2x3 composite system, is separable. On the other hand, 
the general problem of finding sufficient and necessary condition for separability in higher dimensions remains open 
£^ , (see e.g. ||[| and references therein). 

The question of how many mixed quantum states are separable has been raised in |^,||. In particular, it has 
been shown that the relative likelihood of encountering a separable state decreases with the system size TV, while a 
Q/ neighborhood of the maximally mixed state, /?* ~ I/TV, remains separable gj-Q]. 

From the point of view of a possible applications it is not only important to determine, whether a given state is 
entangled, but also to quantify the degree of entanglement. Among several such quantities |§-[lj|], the entanglement of 
formation introduced by Bennet et al. [p"2| is often used for this purpose. Original definition, based on a minimization 
procedure, is not convenient for practical use. However, in recent papers of Hill and Wootters p^JTifl the entanglement 
of formation is explicitly calculated for an arbitrary density matrix of the size TV = 4. 

Any reasonable measure of entanglement have to be invariant with respect to local transformations J^]. In the 
problem of d spin 1/2 particles, for which TV = 2 d , there exist 4 d — 3d+ 1 invariants of local transformations (lq| , and 
' all measures of entanglement can be represented as a function of these quantities. In the simplest case d = 2 there 
exists 9 local invariants, [[l5|-[i"8"|| . These real invariants fix a state up to a finite symmetry group and 9 additional 
discrete invariants (signs) are needed to make the characterization complete. Makhlin has proved that two states are 
locally equivalent if and only if all these 18 invariants are equal Jl^] . Local symmetry properties of pure states of two 
and three qubits where recently analyzed by Carteret and Sudbery j2(3]. A related geometric analysis of the 2x2 
composed system was recently presented by Brody and Hughston pi| . 

The aim of this paper is to characterize the space of the quantum "effectively different" states, i.e. the states 
non equivalent in the sense of local operations. In particular, we are interested in the dimensions and geometrical 
properties of the manifolds of equivalent states. In a sense our paper is complementary to pc| , m which the authors 
consider pure states for three qubits, while we analyze local properties of mixed states of two subsystems of arbitrary 
size. 

We start our analysis defining in section II the Gram matrix corresponding to any density matrix p. We provide 
an explicit technique of computing the dimension of local orbits for any mixed state of the general K x M problem. 
In section III we apply these results to the simplest case of 2 x 2 problem. We describe a stratification of the 6D 
manifold of the pure states and introduce the concept of absolute separability. A list of non generic mixed states of 
TV = 4 leading to submaximal local orbits is provided in the appendix. 
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II. THE GRAM MATRIX 



A. 2 x 2 system 

For pedagogical reasons we shall start our analysis with the simplest case of the 2x2 problem. The local trans- 
formations of density matrices form a six-dimensional subgroup C = SU(2) ® SU(2) of the full unitary group C/(4). 
Let W denote a Hermitian density matrix of size 4, representing a mixed state. Identification of all states which 
can be obtained from a given one W by a conjugation by a matrix from C leads to the definition of the "effectively 
different" states, all effectively equivalent states being the points on the same orbit of SU(2) £g) SU{2) through their 
representative W . 

The manifold W pure of N = 4 pure states, equivalent to the complex projective space, CP 3 , is 6 dimensional. 
Although both the manifold of pure states and the group of local transformations are six-dimensional it does not 
mean that there is only one nontrivial orbit on Wpure- Indeed, at each point W £ W pU re local transformations C/(s), 
parametrized by six real variables s = (si, . . . , sq), such that U(0) equals identity, determine the tangent space to the 
orbit, spanned by six vectors: 

^ : =O s= o = & [7(S)Wt(S)l ^- (1) 

The dimension of the tangent space (equal to the dimension of the orbit) equals the number of the independent Wi 
and, as we shall see, is always smaller then six. Using the unitarity of U(s) one easily obtains 

= {k,W], (2) 

with U :— {dU/dsi) s=0 , and establishes the hermiticity of each Wi. 

Although the so obtained Wi depend on a particular parametrization of C/(s), the linear space spanned by them 
does not. In fact we can choose some standard coordinates in the vicinity of identity for each SU(2) component 
obtaining 

h=i(Jk®I, h+3 = I ®i°k, (3) 

where ct^, k — 1,2,3 stand for the Pauli matrices, and / is the 2x2 identity matrix. Obviously, the antihermitian 
matrices ij, i = 1, . . . , 6, form a basis of the su 2 © su 2 Lie algebra. 

The dimensionality of the tangent space can be probed by the rank of the real symmetric 6x6 Gram matrix 

C mn := h±W m W n . (4) 

formed from the Hilbert-Schmidt scalar products of W^'s in the space of Hermitian matrices. The most important part 
of our reasoning is based on transformation properties of the matrix C along the orbit. In order to investigate them let 
us assume thus, that W' and W are equivalent density matrices, i.e. there exists a local operation U £ SU(2) ® SU (2) 
such that W' = UWU^ . A straightforward calculation shows that the corresponding matrix C calculated at the 
point W' is given by: 

C' mn = ^TrW' m W' n = ^Tr {[l' m , W] [l' n , W]) , (5) 

where 

l'i ■= UUiU, i=l,..., 6. (6) 

The transformation (^) defines a linear change of basis in the Lie algebra su 2 ffi su 2 and as such is given by a 6 x 6 
matrix O i.e. l'i — Ylj=i Oijlj- It can be established that O is a real orthogonal matrix: O -1 = T , cither by the 
direct calculation using some parametrization of SU{2) (g) SU{2) respecting (|^), or by invoking the fact that su 2 ©su 2 
is a real Lie algebra and defines the adjoint representation of SU(2) <g> SU{2). 

Using the above we easily infer that matrices C corresponding to equivalent states are connected by orthogonal 
transformation: C = OCO T . It is thus obvious that properties of states which are not changed under local transfor- 
mations are encoded in the invariants of C, which can thus suit as measures of the local properties such as entanglement 
or distilability. As shown in the following section, the above conclusions remains valid, mutatis mutandis, if we drop 
the condition of the purity of states and go to higher dimensions of the subsystems. 



Wi :-- 



IT) > W 
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B. General case: K x M system 



A density matrix W (and, a fortiori, the corresponding matrix C) of a general bipartite K x M system can be 
conveniently parametrized in terms of (KM) 2 — 1 real numbers dj, b a , Gj a , j — 1, . . . , K 2 — 1, a = 1, . . . , M 2 — 1 as 



W — , „ J + iak(ek ®I) + ib a {I® f a ) + G ka {e k ® f a ), (7) 



1 

where efe and / a are generators of the Lie algebras suk, and sum fulfilling the commutation relations 

[ej,ek] = Cjkiei, [f a , fp] = rf a /3 7 / 7 , (8) 

normalized according to: 

Treje fc = -25 jk , Trf a fp = -25 a p. (9) 

In the above formulas we employed the summation convention concerning repeated Latin and Greek indices. We also 
used the same symbol I for the identity operators in different spaces, as their dimensionality can be read from the 
formulas without ambiguity. Positivity of the matrix W imposes certain constraints on on the parameters Oj , b a , Gj a . 

By analyzing the effect of a local transformation L = V ® U € SU(K) (g> SU(M) upon W we see that a := (etj), 
j = 1, . . . , -fiT 2 — 1 and b = (b a ), a = 1, . . . , M 2 — 1 transform as vectors with respect to the adjoint representations 
of SU(K) and SU(M), respectively, whereas G :— (G^ a ) is a vector with respect to both adjoint representations. 

In analogy with the previously considered case of pure 2x2 states, we can choose the parametrization of the local 
transformations in such a way that the tangent space to the orbit at W is spanned by the vectors 

W i = [e i ®I,W], W a = [I®f a ,W]. (10) 

The number of linearly independent vectors equals the dimensionality of the orbit. As previously this number is 
independent of the chosen parametrization and can be recovered as the rank of the corresponding Gram matrix C, 
which takes now a block form respecting the division into Latin and Greek indices 



C 



A B 
B T D 



(11) 



where 



Aij = ^TrWiWj, B ia = iTrT4W Q , D a0 = ^TvW a Wp. (12) 

The Gram matrix C has dimension K 2 + M 2 — 2, the square matrices A and D are (K 2 — 1) and (M 2 — 1) 
dimensional, respectively, while the rectangular matrix B has size (K 2 — 1) x (M 2 — 1). The matrix C is nonnegative 
definite and the number of its positive eigenvalues gives the dimension of the orbit starting at W and generated by 
local transformations. A direct algebraic calculation gives 

Aij = (2GkaG ma + MakO, m )CiklCj m l: 
^ia 2GkflGy n ~yCik>mdot^fi 1 

D a p = (2G my G m s + Kb^bs)d a7fl dpsn- (13) 
In this way we arrived at the main result of this paper: 

Dimension Di of the orbit generated by local operations acting on a given mixed state W of any K x M bipartite 
system is equal to the rank of the Gram matrix C given by ( pl[ ) - p^j. 

If all eigenvalues of C are strictly positive the local orbit has the maximal dimension equal to Di = K 2 + M 2 — 2. 
In the low dimensional cases it was always possible to find such parameters aj and b a , i.e. such a density matrix W 
that the local orbit through W was indeed of the maximal dimensionality. We do not know if such an orbit exists in 
an arbitrary dimension K x M, although we suspect that is the case in a generic situation (i.e. all eigenvalues of W 
different, nontrivial form of the matrix G). In the simplest case 2 x 2 we provide in Appendix A the list of all, non 
generic density matrices corresponding to sub-maximal local orbits. All other density matrices lead thus to the full 
(six) dimensional local orbits. 

This approach is very general and might be applied for multipartite systems of any dimension. Postponing these 
exciting investigations to a subsequent publication ]22] ], we now come back to the technically most simple case of 
original 2 x 2-dimensional bipartite system. 
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III. LOCAL ORBITS FOR THE 2x2 SYSTEM 



A. Stratification of the 6D space of pure states 



The pure states of a composite 2x2 quantum system form a six-dimensional submanifold W pur e of the fifteen- 
dimensional manifold of all density matrices in the four-dimensional Hilbert space, i.e. the set of all Hermitian, 
non-negative 4x4 matrices with the trace 1. Indeed, the density matrices W and W of two pure states described 
by four-component complex, normalized vectors and |w/)(t//| coincide, provided that \w'} — U\w), where U is 

a unitary 4x4 matrix which commutes with W. Since W has threefold degenerate eigenvalue 0, the set of unitary 
matrices rendering the same density matrix via the conjugation W' = UWU\ can be identified as the six dimensional 
quotient space U(4)/[U(3) x U(l)] = CP 3 . The manifold of the pure states itself is thus given as the set of all 
matrices obtained from Wo := \wo}{wo\, where \wq) = [1,0, 0,0] T , by the conjugation by an element of CP 3 and 
conveniently parametrized by three complex numbers x,y,z: \w) := Af[l, x, y, z] T , W = W(x, y, z) := |w)(u>|, where 
J\f = (1 + |x| 2 + \y\ 2 + Iz] 2 )^ 1 / 2 is the normalization constant, and we allow the parameters to take also infinite values 
of (at most) two of them. In more technical terms we consider thus the orbit of U(4) through the point Wo in the 
space of Hermitian matrices. 

In fact, since the normalization of density matrices does not play a role in the following considerations, we shall 
take care of it at the very end, and parametrize the manifold of pure states by four complex numbers v, x, y, z being 
the components of \w), (the overbar denotes the complex conjugation): 



\w) 



W =\w)(w\ 



vv vx vy vz 

xv xx xy xz 

yv yx yy yz 

zv zx zy zz 



(14) 



bearing in mind, when needed, that the sum of their absolute values equals one. In fact, equating one of the four 
coordinates with a real constant yields one of four complex analytic maps which together cover the complex projective 
space CP 3 (with which the manifold of the pure states can be identified) via standard homogeneous coordinates. This 
leads to a more flexible, symmetric notation, and dispose off the need for infinite values of parameters. 

The dimensionality of the orbit given by rank(C) is the most obvious geometric invariant of the orthogonal trans- 
formations of C. As it should it does not change along the orbit. All invariant functions (or separability measures) 
can be obtain in terms of the functionally independent invariants of of the real symmetric matrix C under the action 
of the adjoint representation of SU(2) ® SU(2). In particular, the eigenvalues of C are, obviously, such invariants. 
Substituting our parametrization of pure states density matrices ( |l4]) to the definition of C (^) yields, after some 
straightforward algebra, the eigenvalues 



Ai = 0, A 2 



A 3 = A 4 = 1 + 2M, A 5 = A 6 = 1 - 2\w\ 



(15) 



where lo := vz — xy. For any pure state one may explicitly calculate the entropy of entanglement fl2|| or a related 
quantity, called concurrence jjj. For the pure state (^4|) the concurrence equals 



c = 2\uj\ = 2\vz — xy\ 
and c £ [0, 1]. Thus the spectrum of the Gram matrix may be rewritten as 

eig(C) = {0, 2c 2 , 1 + c, 1 + c, 1 - c, 1 - c}. 



(16) 



(17) 



The number of positive eigenvalues of C determines the dimension of the orbit generated by local transformation. 
As already advertised, the dimensionality of the orbit is always smaller then 6. In a generic case it equals 5, but for 
lo = (c = - separable states) it shrinks to 4 and for |o;| = 1/2 (c = 1 - maximally entangled states) it shrinks to 3. 
These results have already been obtained in a recent paper by Carteret and Sudbery [ p0[ , who have shown that the 
exceptional states (with local orbits of a non-generic dimension) are characterized by maximal (or minimal) degree of 
entanglement. 

In order to investigate more closely the geometry of various orbits let us introduce the following definition: 



Wn ■■= {W = ww f : w = [v, x, y, zf e C 4 , ||w|| 2 = \v\ 
+ |y| 2 + |z| 2 = l, \(vz-xy)\ = n}. 



(18) 



4 



It is also convenient to define a map from the space of state vectors {w = [v,x,y, z] T £ C 4 : ||w|| 2 = \v\ 2 + \x\ 2 + 
\y\ 2 + \z\ 2 = 1} to the space of complex 2x2 matrices 



X(w) 



v y 

x z 



(19) 



In terms of X(w) the length of a vector w and the bilinear form w(w) 
and w(w) = detX(w). From the Hadmard inequality 



vz - 



|detX(w)|<[(|t;| 2 + N 2 )(M 2 + |z| 2 )] 1 /2 ; 



■ xy read thus: ||w|| 2 = TrX(w)Xt(w) 



(20) 



we infer |w(w)| < ^. Indeed, since |t>| + \x\ + \y\ + \z\ = 1, the right-hand-side of (|20| ) equals its maximal value of 
2 for |t>| 2 + \x\ 2 = | = \y\ 2 + \z\ 2 . A straightforward calculation shows also, that a local transformation L = V <g> U 
sends w to w' = Lw if and only if X(w') = UX(w)V T . As an immediate consequence we obtain the conservation of 
|oj(w)| under local transformation. Together with the obvious conservation of | |w|| (which, by the way, is also easily 
recovered from ||w|| 2 = TrX (w)X^ (w)), it shows that the parametrization ( |l8| ) is properly chosen. Moreover it can 
be proved that C acts transitively on submanifolds ( |Tg| ) of constant \u>\, i.e. for each pair W — ww', W — w'w'^ 
such that |w(w)| = |lj(w')| = f2, there exists such a local transformation L 6 C that W 1 = L(W) := LWL\ or, 
in other words, that the manifold (|l8|) of constant |w| is an orbit of the group of local transformations £ through a 
single point W i.e. Wn — C{W). To this end, it is enough to show that each W € Wn can be transformed by a local 
transformation into Wg = wgwj, where wg = [cos(6'/2), 0, 0, sin((9/2)] T with sin9 — 2ui (from the above mentioned 
bound for |w(w)| we know that it is sufficient to consider < 9 < ^). To this end we invoke the singular value 
decomposition theorem which states that for an arbitrary (in our case 2x2) matrix X, there exist unitary U',V 
such that 



X 1 := U'XV' 



p 
q 



p>q>0. 



Let now V = e^V , U' = e ir >U, V, U S SU(2). We can rewrite (||) as 



UXV T = 



pe 




i4 Q 



qe 



i(f> 



p>q>0cf>:=-(r) + 



(21) 



(22) 



< < \ connects 



Substituting X = X(w) (|19|), we obtain p 2 + q 2 = TyXX^ = ||w|| 2 = 1 and invoking the invariance of pq = |detX| = 
|oj(w)| = sin 29. This gives an unique solution p = cos 9, q = sin9 in the interval < 9 < 5. On the other hand, as 
above mentioned, the transformation ( p2| ) corresponds to Iw = w' g = [cos(6*/2)e l< ^, 0, 0, sin(#/2)e^] T , but obviously 
Wg — w e w ^ = Wg = wjwj, i.e., finally, LWL^ — Wg, with L = V <£> U e C as claimed. This is, obviously, a 
restatement of the Schmidt decomposition theorem for 2 x 2 systems. 

Now we can give the full description of the geometry of the states. The line into Wg = wgw 
all "essentially different" states. At each 9 different from 0,7r/2 it crosses a five-dimensional manifold of the states 
equivalent under local transformations. The orbits of submaximal dimensionality correspond to both ends of the 
line. For 9 — ir/2 the orbit is three-dimensional. The states belonging to these orbits are maximally entangled, since 
\oj\ — 1/2 corresponds to c = 1. 

In order to recover the whole orbits we should find the actions of all elements of the group of local transformations 
on a representative of each orbit (e.g. one on the above described line). Since, however, the orbits have dimensions 
always lower than the dimensionality of the group, the action is not effective, i.e. for each point on the orbit, there is 
a subgroup of C which leaves this point unmoved. This stability subgroups are easy to identify in each case. Taking 
this into account we end up with the following parametrization of three-dimensional orbits of the maximally entangled 
states 



W V4 = {W 



w(a,xi,X2)}, w(a,xi,x 2 )} = -j= 



cosae 
sin ae K 



sin ae 



->X2 



cos ae 



-i-Xl 



(23) 



with < Xi < 2?r, < a < tt/2, which means that topologically this manifold is a real projective space, MP 3 = S* 3 /^, 
where Z2 is a two elements discrete group. This is related to the well known result that for bipartite systems the 
maximally entangled states may be produced by an appropriate operation performed locally, on one subsystem only. 
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The manifold of maximally entangled states (|23|) is cut by the line of essentially different states at the origin of the 
coordinate system (a,xi,X2)- 

The four dimensional orbit corresponding to 9 — consists of separable states characterized by the vanishing 
concurrence, c = 0. The parametrization of the whole orbit, exhibiting its S 2 x S 2 structure, is given by: 



w(a,/3,Xi,X2) 



cos a cos [3e tXl 
cos a sin (3e lX2 
sin a cos /3e~' iX2 
sinasin/3e~ iXl 

< Xt < 2tt, < a,P < 7r/2. 



(24) 



The majority of states, namely these which are neither separable, nor maximally entangled, belongs to various 
five-dimensional orbits labeled by the values of the parameter 9 with < 9 < ir /2. In this way we have performed a 
stratification of the 6D manifold of the pure states, depicted schematically in Fig. lb. 

For comparison we show in Fig. la the stratification of a sphere S 2 , which consists of a family of ID parallels and 
two poles. Zero dimensional north pole on CP 1 corresponds to the 3D manifold of maximally entangled states in 
CP 3 , while the 4D space of separable states may be associated with the opposite pole. In the case of the sphere (the 
earth) the symmetry is broken by distinguishing the rotation axis pointing both poles. In the case of N = 4 pure 
states the symmetry is broken by distinguishing the two subsystems, which determines both manifolds of maximally 
entangled and separable states. 




b)N = 4 



(maximally 
entangled) 




oo > 



CD = 

(separable) 



V = CV 



FIG. 1. Stratification of the sphere along the Greenwich meridian (a), stratification of the 6 dimensional space of the N = 4 
pure states along the line of effectively different states, ui £ [0, 1/2] (b). The poles correspond to the distinguished submanifolds 
of CP 3 : the 3D manifold of maximally entangled states and the 4D manifold of separable states. 



G 



B. Dimensionality of global orbits 



Before we use the above results to analyze the dimensions of local orbits for the mixed states of the 2x2 problem, 
let us make some remarks on the dimensionality of the global orbits. The action of the entire unitary group U{4) 
depends on the degeneracy of the spectrum of a mixed state W. Let W — VRV^ , where V is unitary and the diagonal 
matrix R contains non negative eigenvalues r£. 

Due to the normalization condition TrW = 1 the eigenvalues satisfy r% + T% + r% + = 1. The space of all 
possible spectra forms thus a regular tetrahedron, depicted in Fig. 2. Without loss of generality we may assume that 
?*i > T2 > T"3 > f4 > 0. This corresponds to dividing the 3D simplex into 24 equal asymmetric parts and to picking one 
of them. This set, sometimes called the Weyl chamber ^3|, enables us to parametrize entire space of mixed quantum 
states by global orbits generated by each of its points. 



a) (0001) 




FIG. 2. The simplex of eigenvalues of the N = 4 density matrices (a). Pure states are represented by four corners of the 
thetrahedron, while its center denotes the maximally mixed state p». Magnification of the asymmetric part od the simplex, 
related to the Weyl chamber (b). It can be decomposed into 8 parts according to different kinds of degeneracy of the spectrum. 

Note that the unitary matrix of eigenvectors V is not determined uniquely, since W = VRV^ = VHRWV* , where 
H is an arbitrary diagonal unitary matrix. This stability group of U is parametrized by N — 4 independent phases. 
Thus for a generic case of all eigenvalues r» different, (which corresponds to the interior Kim of the simplex), the 
space of global orbits has a structure of the quotient group C/(4)/[[/(l) 4 ]. It has D g = 16 — 4 = 12 dimensions. 
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If degeneracy in the spectrum of W occurs, say r\ — r 2 > r 3 > r 4 , than the stability group H = U(2) x U(l) x U(l) 
is 4 + 1 + 1 = 6 dimensional p4|. In this case, corresponding to the face K211 of the simplex, the global orbit U/H 
has D g = 16 — 6 = 10 dimensions. The dimensionality is the same for the other faces of the simplex, K\n and 
Kxxi- The important case of pure states corresponds to the triple degeneracy, t\ > r 2 = r 3 = r 4 for which the 
stability group H equals Z7(3) x U(l). The orbits U/H = SU(4)/U(2>) have a structure of complex projective space 
CP 3 . This 6D manifold results thus of all points of the Weyl chamber located at the edge K13. These parts of the 
asymmetric simplex are shown in Fig. 2, the indices labeling each part give the number of degenerated eigenvalues in 
decreasing order. For another edge K22 of the simplex H = U(2) X [7(2) and the quotient group U/H is 16 — 8 = 8 
dimensional. In the last case of quadruple degeneracy, corresponding to the maximally mixed state = 1/4, the 
stability group H = Z7(4), thus Di = 0. A detailed description of the decomposition of the Weyl chamber with respect 
to the dimensionality of global orbits for arbitrary dimensions is provided in p5fl. 



C. Dimensionality of local orbits 



For K = M = 2 (two qubit system) c^fc 
tensor. Formulae (113) give in this case 



2e ijk and 



Pi 



-2e Q| g 7 , where e a( a 7 is completely antisymetric 



A = 8[(TrG'G' J ) • I - G G ] + 8(||a' 
D = 8[(TrG'G' T ) • / - G' T G'} + 8(||b 



,'l|2 



a'a' T ) 



JU2 



WW 



and 



BG' T = G' T B = -16 detG' ■ I. 



(25) 
(26) 

(27) 



where 3D vectors a', b' and a 3 x 3 matrix G' represent a certain N = 4 mixed state W in the form (0). For 
later convenience we denote the system variables by symbols with primes. For detG' =/= the last equation gives 
B = — 16detG' T • (G but below we will show the more convenient representation of B. 

Since G' is real, we can find its singular value decomposition in terms of two real orthogonal matrices 0%, O2 and 
a positive diagonal matrix 



OxG'Ol = G = 



m 

(12 
^3 



Mi > A*2 > M3 > 0. 



(28) 



If the determinant of G' is positive then one can choose 0\ and O2 as proper orthogonal matrices (i.e. with the 
determinants equal to one). In this case the singular value decomposition (|8|) corresponds to a local transformation 
W = U\ ® U2W(U\ ® U2) 1 ' ■ In the opposite case of a negative determinant of G' one of the matrices 0\, O2 has also 
a negative determinant. Alternatively, we can assume that 0\ and O2 are proper orthogonal matrices (with positive 
determinants), and, consequently, the singular value decomposition corresponds to a local transformation, but with 
Mi < A*2 < M3 < 0. 

From ( p6| ) and (|2^) it follows, that the above transformation G = OiG'O^j if supplemented by a = Oia' and 
b = 2 b', induces the transformation G = G'(G', a', b') ^ G(G, a, b) = {Ox 2 )G'(Oi © 2 ) T , where 



Oi © 2 := 



Oi 
2 



leaving the spectrum of G invariant. The explicit form of the transformed matrix inferred from (26) 
reads 



(29) 

ID, and (H) 



G = 



8 (f4 + a*§) 



+16/x 2 /i 3 









Mi) 



8(M? 



+ 16^iipi3 





8(M? 



Ml) 



+ 16/x 2 /z 3 





(mI 



m1) 



8(||a| 



| 2 -/-a(a) T ) 




+ 16/xi/x 2 


8(||b|| 2 J-b(b) 





+16/ii/i 3 



8 (/. 2 + 







+ 16/xi/i 2 



8 ( M 2 + ^ 2 ) 



:=G 



G 



G 



a.b 



(30) 
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which is the sum of two real positive definite matrices, Cq and C a .b- Their eigenvalues are, respectively 

pi = 8(/ii + /i 2 ) 2 , p 2 = 8(mi + M3) 2 , P3 = 8(A*2 + M3) 2 , 

Pa = 8(/ii - /i 2 ) 2 , p 5 = 8(/ii - M3) 2 , Pe = 8(/x 2 - M3) 2 , (31) 

and 

t>i = v 2 = ||a|j 2 , v 3 = z/ 4 = ||b|| 2 , z/ 5 = ^ 6 = 0. (32) 

Although two parts, Cg and C ai b of C, usually, do not commute and the eigenvalues Ai > ■ • • > Xq > of C cannot 
be immediately found, we can investigate the possible orbits of submaximal dimensionalities using the fact that both 
Cq and C a b are positive definite. It follows thus that the number of zero values among the eigenvalues Ai, . . . , Xq 
of C has to be matched by at least the same number of zeros among p\, . . . ,p$ and among t/i,...,i>e, moreover the 
eigenvectors to the zero eigenvalues of the whole matrix C are also the eigenvectors of the components Cg and C aj b 
(also, obvoiusly, corresponding to the vanishing eigenvalues) 

The co-rank r' c G (the number of vanishing eigenvalues) of Cg equals 

6 for pi = p 2 = M3 — G = 0, 

3 for pi = P2 = P3 '■= M 7^ & G = pi, 

2 for p:= pi> p 2 = P3 = 0, 

1 for p M ■= Pi > p-2 = P3 ■= Mm ^ 0, 

or ^ p := pi = p 2 > P3, (33) 

and is equal in all other cases, whereas for C a .b it co-rank r'c a b reads 

6 for a = b = 0, 

4 for a = 0, b ^ or a ^ 0, b = 0, 

2 for a ^ 0, b ^ 0. (34) 

As already mentioned, in a generic case all eigenvalues of the 6D Gram matrix C are positive and the dimension of 
local orbits is maximal, di — 6. On the other hand, the above decomposition of the Gram matrix is very convenient to 
analyze several special cases, for which some eigenvalues of C reduce to zero and the local orbits are less dimensional. 
To find all of them one needs to consider 9 combinations of different ranks of the matrices Cg and C aj b as shown in 
the Appendix. 

For any point of the Weyl chamber we know thus the dimension D g of the corresponding global orbit. Using above 
results for any of the globally equivalent states W (with the same spectrum) we may find the dimension Di of the 
corresponding local orbit. This dimension may be state dependent, as explicitly shown for the case of TV = 4 pure 
states. Let D m denotes the maximal dimension Di 7 where the maximum is taken over all states of the global orbit. The 
set of effectively different states, which cannot be linked by local transformations has thus dimension D e i = D g — D m . 
For example, the effectively different space of the N = 4 pure states is one dimensional, Dd — 6 — 5 = 1. 



D. Special case: triple degeneracy and generalized Werner states 

Consider the longest edge, -RT13, of the Weyl chamber, which represents a class of states with the triple degeneracy. 
They may be written in the form p x := x\^)(^\ + (1 — x)p*, where \^>) stands for any pure state and x £ (0, 1]. The 
global orbits have the structure U(4)/[U(3) x U(l)}, just as for the pure states, which are generated by the corner of 
the simplex, represented by x = 1. Also the topology of the local orbits do not depends on x, and the stratification 
found for pure states holds for each 6 dimensional global orbit generated by any single point of the edge. 

Schematic drawing shown in Fig. 1 is still valid, but now the term " maximally entangled" denotes the entanglement 
maximal on the given global orbit. It decreases with x as for Werner states, with chosen as the maximally 
entangled pure state p7j . For these states the concurrence decreases linearly, c(x) — (3x — l)/2 for x > 1/3 and is 
equal to zero for x < 1/3. Thus for sufficiently small x (sufficiently large degree of mixing) all states are separable, 
also these belonging to one of the both 3D local orbits. This is consistent with the results of ||, where it was proved 
that if Trp 2 < 1/3 the 2x2 mixed state p is separable. 
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This condition has an appealing geometric interpretation: on one hand it represents the maximal 3D ball inscribed 
in the tetrahedron of eigenvalues, as shown in Fig. 3. On the the other, it represents the maximal 15D ball _Bj\/ , (in 
sense of the Hilbert-Schmidt metric, D 2 HS {px, p 2 ) = Tr(pi — P2) 2 ), contained in the 15D set of all mixed states for 
N = 4. Both balls are centered at the maximally mixed state p* (the center of the eigenvalues simplex of side y/2), 
and have the same radius l/2\/3- A similar geometric discussion of the properties of the set of 2 x 2 separable mixed 
states was recently given in [|6| . 



N=4 (0001) 




(0100) 

FIG. 3. Separability of the maximal 15D ball: all mixed states with spectra represented by points inside the ball inscribed 
in the 3D simplex of eigenvalues of the N = 4 density matrices are separable. 

To clarify the structure of effectively different states in this case we consider generalized Werner states 

p(x,0):=x\*o){*o\ + 0--*)P*> ( 35 ) 

where the state \tyg) := [cos(0/2), 0, 0, sin(0/2)], contains the line of effectively different pure states for 9 £ [0, 7r/2]. 
Note that the case 9 = tt/2 is equivalent to the original Werner states |27J]. Entanglement of formation E for the 
states p(x,9) may be computed analytically with help of concurrence and the Wootters formula |Q. The results 
are too lengthy to be reproduced here, so in Fig. 4 we present the plot E = E(x,9). The graph is done in polar 
coordinates, so the pure states are located at the circle x — 1. For each fixed x the space of effectively different states 
is represented by a quarter of the circle. For x < 1/3 entire circle is located inside the maximal ball Bm, and all 
effectively different states are separable. Points located along a circle centerd at p* represent mixed states, which are 
described by the same spectrum and can be connected by a global unitary transformation J7(4). In accordance to the 
recent results of Hiroshima and Ishizaka pq , the original Werner states enjoy the largest entanglement accesible by 
unitary operations. 

The convex set S of separable states contains a great section of the maximal ball and touches the set of pure states 
in two points only. The actual shape of S (at this cross-section) looks remarkably similar to the schematic drawing 
which appeared in Moreover, the contour lines of constant E elucidate important feature of any measure of 
entanglement: the larger shortest distance to <S, the larger entanglement Even though we are not going to prove 
that for any state p, its shortest distance to S at the picture is strictly the shortest in the entire 15D space of mixed 
states, the geometric structure of the function E = E(x,9) is in some sense peculiar: The contours E =const are 
foliated along the boundary of S, while both maximally entangled states are located as far from S, as possible. 
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x E 




-1 -0.5 0.5 1 



FIG. 4. Entanglement of formation E for the generalized Werner states p Xt e represented in the polar coordinates. Intersection 
with the maximal ball centered at p* is separable (white). Dashed horiznotal line, joining two maximally entangled states (*), 
(black), represents the original Werner states. Entanglement E of a mixed state p may be interpreted as its distance from the 
set of separable states. 

E. Absolutely separable states 

Defining separability of a given mixed state p, we implicitly assume that the product structure of the composite 
Hilbert space is given, Ti. — Ti.^ <X> TLb- This assumption is well justified from the physical point of view. For example, 
the EPR scenario distinguishes both subsystems in a natural way ('left photon' and 'right photon'). Then we speak 
about separable (entangled) states, with respect to this particular decomposition of H. Note that any separable pure 
state may be considered entangled, if analyzed with respect to another decomposition of Ti.. 

On the other hand, one may pose a complementary question, interesting merely from the mathematical point of 
view, which states are separable with respect to any possible decomposition of the N = K x M dimensional Hilbert 
space TL. More formally, we propose the following 

definition. Mixed quantum state p is called absolutely separable, if all globally similar states p' = UpU^ are 
separable. 

Unitary matrix U of size N represents a global operation equivalent to a different choice of both subsystems. It is 
easy to see that the most mixed state is absolutely separable. Moreover, the entire maximal ball Bm — B{p*, l/2\/3) 
is absolutely separable for N = 4. This is indeed the case, since the proof of separability of Bm provided in ^ relays 
only on properties of the spectrum of p, invariant with respect to global operations U. Another much simpler proof 
of separability of Bm follows directly from inequality (9.21) of the book of Mehta J29|. 

Are there any 2x2 absolutely separable states not belonging to the maximal ball Bm? Recent results of Ishizaka 
and Hiroshima |30) suggest, that this might be the case. They conjectured that the maximal concurrence on the local 
orbit determined by the spectrum {n, r2, ^3, r^} is equal to c* = max{0, r% — — l^fr^xX^. This conjecture has been 
proved for the density matrices of rank 1,2 and 3 po| . If it is true in the general case than the condition c* > 
defines the 3D set of spectra of absolutely separable states. This set belongs to the regular tetrahedron of eigenvalues 
and contains the maximal ball Bm- For example, a state with the spectrum {0.47,0.30,0.13,0.10} does not belong 
to Bm but its c* is equal to zero. 

IV. CONCLUDING REMARKS 

In order to analyze geometric features of quantum entanglement we studied the properties of orbits generated by 
local transformations. Their shape and dimensionality is not universal, but depends on the initial state. For each 
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quantum state of arbitrary K x M problem we defined the Gram matrix G, the spectrum of which remains invariant 
under local transformation. The rank of G determines dimensionality of the local orbit. For generic mixed states the 
rank is maximal and equal to Di = K 2 + M 2 — 2, while the space of all globally equivalent states (with the same 
spectrum) is {KM) 2 — KM dimensional. Thus the set of states effectively different, which cannot be related by any 
local transformation, has D d = {KM) 2 - KM - {K 2 + M 2 - 2) dimensions. 

For the pure states of the simplest 2x2 problem we have shown that the set of effectively different states is one 
dimensional. This curve may be parametrized by an angle emerging in the Schmidt decomposition: it starts at a 3D 
set of maximally entangled states, crosses the 5D spaces of states of gradually decreasing entanglement, and ends at 
the 4D manifold of separable states. 

We presented an explicit parametrization of these submaximal manifolds. Moreover, we have proved that any pure 
state can be transformed by means of local transformations into one of the states at this line. In such a way we 
found a stratification of the 6D manifold CP 3 along the line of effectively different states into subspaces of different 
dimensionality. 

Since for N = 4 pure states the set of effectively different states is one dimensional, all measures of entanglement 
must be equivalent (and be functions of, say, concurrence or entropy of formation). This is not the case for generic 
mixed states, for which Dd = 6. Hence there exist mixed states of the same entanglement of formation with the same 
spectrum (globally equivalent), which cannot be connected by means of local transformations. 

It is known that some measures of entangled do not coincide (e.g. entanglement of formation E and distillable 
entanglement Ed 11 1). To characterize the entanglement of such mixed states one might, in principle, use 6 suitably 
selected local invariants. This seem not to be very practical, but especially for higher systems, for which the dimension 
Dd, of effectively different states is large and the bound entangled states exist (with E d = and E > 0), one may 
consider using some additional measures of entanglement. All such measures of entanglement have to be functions of 
eigenvalues of the Gram matrix G or other invariants of local transformations |l(]|l^,[l7] |l9| . 

We analyzed geometry of the convex set of separable states. For the simplest N — 4 problem it contains the 
maximal 15D ball, inscribed in the set of the mixed states. It corresponds to the 3D ball of radius 1/2 V3 inscribed 
in the simplex of eigenvalues. This property holds also for 2 x 3 problem, for which the radius is l/\/30. For larger 
problems K X M = N > 8, it is known that all mixed states in the maximal ball (of radius {N{N — l)) -1 ' 2 ) are not 
distillable but the question whether they are separable remains open. 
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APPENDIX A: SUBMAXIMAL LOCAL ORBITS FOR 2x2 PROBLEM 



In this appendix we give the list of all possible submaximal ranks of the Gram matrix G which determine the 
dimension of the local orbit D\ — 6 — Tq- The symbol r' x denotes the co-rank, it is the number of zeros in the 
spectrum of X. In each submaximal case we provide the density matrix W , Gram matrix G and its eigenvalues A,, 
i = 1, . . . , 6 expressed as a function of the the singular values of the matrix G" and the vectors a = Oia' and b = 02 b', 
where orthogonal matrices 0\ and O2 are determined by the singular value decomposition of G . 

In a general case the density matrix W = W (G, a, b) = W {fii,fX2, A*3, ffli, 0-2, 03, b\, b 2 , ^3) is given by 



W 



-a 3 - 63 - jU 3 -61 - ib 2 -ai - ia 2 -/xi + /i 2 

-b x + ib 2 -a 3 + b 3 + ^ 3 -/ix - fi 2 -ai - ia 2 

-a! + ia 2 - [i 2 a 3 - 63 + /i 3 -6x - ib 2 

-Mi + —a-i + ia 2 — 6j + ib 2 a 3 + b 3 — fi 3 



(AI) 



where we use the rotated basis in which G is diagonal. The characteristic equation of the density matrix W reads 



det {W - g) 



Q 4 - Q 3 + 



(» 



|-2 W 



2||b|| -2TrG 



2TrG 4 - (TrG 2 ) 



Q 

2 1 , 



1 

16 



a|| 2 4l|b|| 2 



a|r + ||b| 
1, 



TrG 2 + 8aGb - 8 det G g 



-TrG -2aGb + 2detG 
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-4 ||Ga|| 2 - 4 ||Gb|| 2 + 2 (||a|| 2 + ||b|| 2 ) TrG 2 + 8 (axh^s + a 2 b 2 ^3 + a 3 b 3 ^ 2 ) + (A2) 

It is interesting to note that the characteristic equation of the partially transposed matrix W — W T2 differs only by 
signs of three terms: 



det (W - q) = - g 3 



jj-2||a|| 2 -2||b|| 2 -2TrG 2 



-2 

Q 



1 

16 



TrG 2 + 8aGb + 8 det G )g 



Ia|| 2 -||b|| 2 



2TrG 4 - (TrG 2 ) - - llall 2 - - llbll 2 - -TrG 2 - 2aGb - 2dctG 
V j g 11 11 g 11 11 8 



4 \\Gaf - 4 ||Gb|| z + 2 (J|af + ||b|| z J TrG - 8 (aihtolMi + a 2 b 2 ^3 + a 3 b 3 ^ lf i 2 ) + — . (A3) 



Let Qi and Qi,i= 1,2, 3, 4, denote the eigenvalues of W and W, respectively. Due to Peres-Horodeccy partial transpose 
criterion |fl|p| positivity of Qi may be used to find, under which conditions W is separable. 

In order to compute the concurrence of the density matrix W, let us define an auxiliary hermitian matrix 



W := Wa 2 <g> a 2 W*a 2 <g> a 2 , 



(A4) 



where * represents the complex conjugation. Let £f, i = 1, 2, 3, 4 denote the eigenvalues of W, arranged in decreasing 
order. Then the concurrence c of W is given by [ p~3| , p~4[ 

(A5) 



:= max ^0, - a/6 - V6 - \f£ 



The Gram matrix G = G (G, a, b) = C (fix, /12, A«3, ai, 02, 03, b%, b 2 , 63) corresponding to the density matrix W , reads 
in the general case 



G = 8 



—axa 2 
-axa 3 
-2/i2A*3 






— aia 2 
a? + o§ + Mi + Ms 
-a 2 a 3 


-2/il/i3 





-axa 3 
—a 2 a 3 
a\+al+ /if + fi% 





-2/i X /i 2 



-2/i 2 M3 

-2M1M3 

-2mn2 

b 2 + b i + ^ 2 +^l -b\b 2 -bxh 

-&1&2 b 2 1 + bj + fi 2 1 + 14 -b 2 b 3 

-6163 -b 2 b 3 bl + b 2 2 +fij+ & 

(A6) 



Below we provide a list of the classes of states corresponding to the submaximal ranks tq of the Gram matrices. 
The list is ordered according to the increasing dimensionality of local orbits; Di — rc = 6 — r' c . 

Case 1. r' c = 6, G = 0, a = 0, b = 0,; G = 



•^1,2,3,4,5,6 — 0; W — —I] 01,2,3.4 — - — 01,2,3,4, £l,2,3,4 — — , 

4 4 16 



thus W is separable and concurrence, c, is equal to zero. 
Case 2. r' c = 4, G = 0, a ^ 0, b = 0, 

Ai,2 = 8||a|| 2 , A3, 4, 5, 6 



(A7) 



(A8) 



£-i,2 = j + ||a||,03 )4 = j- ||a|| 



1 „ , _ 1 „ „ 
gi,2 = ^ + ||a|| ,g 3 ,4 = 1 - ||a|| 



1 2 
6,2,3,4 = T77 a , thus c = 0. 
16 



(A9) 
(A10) 
(AH) 
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W represents a density matrix for ||a|| < \ and then is separable (W > 0). 
Case 3. r' c = 3; G = [il, a = 0, b = 0, 

Ai, 2 ,3 = 32 M 2 , A 4 , 5 ,6 = (A12) 

£•1,2,3 = j - M, = j + 3,u, (A13) 

1 ~ 1 /. X 

£•1,2,3 = -+n, Qi = -- i^. (A14) 

a = ^(12 M +l) 2 , 6,3,4 = ^ (1 - 4 M ) 2 , (A15) 



(A16) 



for ^ < 
^ ~ 2" for T5 - ^ < | 

14/ > for < (i < j and W is separable for |^| < ^. 

Case 4. = 2. G = 0, 

Ai, 2 = 8||a|| 2 , A 3 , 4 - 8||b|| 2 , A 5 , 6 = 0. (A17) 
Qi = J + ||a|| + ||b||, e2 = l-||a||-||b|| ) ft3 = i + |||a||-||b|||, Qi = ± - |||a|| - ||b||| (A18) 

ei = J + ||a|| + ||b||, 2 = i-|N|-||b||, e3 = ^ + |||a||-||b|||, ^4 = J-|||a||-||b|||. (A19) 

a,2 = ^ + (||a|| + ||b!|) 2 , 6 ,4 = ^ + (||a||-||b||) 2 ; c = . (A20 ) 

> for ||a|| + ||b|| < j and is then separable. 
Case 5. r' c = 2 

G = diag(jj,, 0,0), a=[a,0,0] T , b=[fo,0,0] T , 

Ai,2 = 8(a 2 + M 2 ), A 3 , 4 - 8 (b 2 + M 2 ) , A 5 , 6 = (A21) 

£>i = i + a + fe-^, g 2 = 2- a + b + fi, Q 3 = ^- a- b-/j,, p 4 = i+ a- & + ^, (A22) 
g!^2+a + b-fi, g 2 = -^-a + b + LL, g 3 = ^- a- b~iJ,, p 4 = i + a - 6 + /i. (A23) 

6,2= Q+m) -(a-&) 2 , 6,4= (J-m) -(a + fc) 2 , c = 0. (A24) 

> for a = ||a|| < \, b = ||b|| < \, \fi\<\; then W is separable. 

Case 6. r'c = 1 

G = diag(fi, 0, 0), a=[o, 0, 0] T , 
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A 2 = 4 (j|b|| 2 + /i 2 - ^-||b|| 2 ) 2 +V&?) 
A 3 = 8 (||b|| 2 + fi 2 ) , A 4 , 5 = 8 (a 2 + ^ 2 ) , A 6 = 



Qi = j + a + yV + ll b l| 2 - 2& iM> Q2 = 2+a- \jp? + ||b|| 2 - 26^, 
P3 = i-a+ V / A i 2 + ||b|| 2 + 26 lM , W = J-a- V / M 2 + l|b|| 2 + 26 lA i, 



1 



& = - - a + Vm 2 + ||b|r + 26i/i, Qi = j - a - \j [i 2 + ||b|| 2 + 26 lM , 



6,2 - ^ + M 2 - a 2 - ||b|| 2 + V /4( a 2- M 2)||b|| 2 + V&2 + 2M06!, 
6,4 = ^ + M 2 - a 2 - ||b|| 2 - V / 4( a 2 -M 2 )||b|| 2 + 4 / i262 + 2 Ma 6 1 , 
so c — 0. If W represents a density matrix (W > 0) then it is separable. 



Case 7. r' c = 1 

G = fil, b = £a, 



Ai = 4 ((e 1) ||aj| 2 + yjf + ^16^ + I 1} 



2 "aP 



A 2 = 4 ^(C 2 - 1) ||a|| 2 + _ ^ 16m 4 + ^2 _ ^ || a ||2j ; 
As = 4 (V - 1) ||a|| 2 - \Jn 2 + ^lQn A + (.e-l?\\B\A , 



A4 = 4(^ 2 -l)||a|| 2 -V^ 2 
A 5 = 32^ 2 ,A 6 = 



yW + ^-ifllall 2 ^ , 



ei = \ ~ M + |£ + 1| ||a|| , Q2 = \-H-\Z + 1| l|a|| , 



,, = 1 + + V /4 M 2 + (e _i) 2 ||a|| 2 , Q4 = I + M - ^V + (£-l) 2 ||a|| 2 



ei = i+ M +| € -l|||a||, ^2 = \ + M - |£ - 1| l|a|| , 

& = J - M + vV + (£ + i) 2 ||a|| 2 , ^ = J - n - \/v + M 



2 n n 2 

a 



a = ^ + f + 5^ 2 - K - 1) 2 Nl 2 + mVW + i) 2 -i6(£-i) 2 llai 



6 = 1 + I + 5 M 2 - (£ - I) 2 l|a|| 2 - M V 4 (// I V 10 (i Vf ||a|| 



1 



2 



2 n n 2 



£s,4= (4-M -((£ + in|a 
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If W > i.e. Qi > 0, i — 1,2,3,4 then |/x| < | and Qi > 0, i — 1,3, hence T4 7 is nonseparable for 
yV + ^ + lfNI 2 > i - M > IC + 1| l|a|| or I < |£ - l|||a||. 
Case 8. r' c = 1 

G = diag((n,(i 2 ,to), a=[a, 0, 0] T , b=[6, 0,0] T , 



Ai, 2 = 4 (V + 6 2 + + 2^ 2 + y / 16 M 2 M 2 + (a 2_ 6 2 ) 2^ 

.2 i 1,2 , „2 , o,,2 _ ,/ifi„2„2 



A 3 , 4 - 4 (^a 2 + & 2 + 2/4 + 2/4, - ^ 16^?^l + (a 2 - & 2 ) J , A 5 = 32/i 2 , A 6 = (A33) 
2i = | - Mi + a + b , Q2 = j - Hi - a - b, 

Qz = \+ (ii + \j^l + {a-b)\ Qi = - A +(ii- \J^l + (a-bf (A34) 
Qi = i + (jli - a + b, Q 2 = ^+ Hi + a - b, 



i r~, — ; — ^ _ i 



03 = i ~ Mi + V 4 mI + (a + b) , Qi = - A -(ii- \/A(il + (a + b) (A35) 



[\ /':) -{a + bf, &=(t/Q+Mi) - (a - &) 2 + 2 M 2 



J+Mil -(a-&) 2 -2 Al2 | , (A36) 

If > i.e. £>j > 0, i = 1, 2, 3, 4 then |/xi | < \ and & > 0, i = 1, 3, hence W is nonseparable for 4/i| + (a + 6) 2 > 
i - ,ui > |a + b\ or | < 6 - a - (i\. 
Case 9. r' c = 1 



G = diag((ii,(ii,(i 2 ), a =[0,0, a] T , b=[0, 0,6] T . 

Ai, 2 = 4 (V + 6 2 + 2(ij + 2 M 2 + ^/l6/i 2 /^ + (a 2 -6 2 ) 2 ^ , 



A 3 , 4 = 4 (^a 2 + V + 2 M 2 + 2 M 2 - ^ 16/i 2 /; 2 + (a 2 - 6 2 ) z ) , A 5 = 32 M f, A 6 = 0, (A37) 
Qi = ^ - M2 + a + &, £>2 = j - M2 - a - b, 



03 



\+to + yjtfil + ia-b) 2 , Qi = - A +(i 2 -\J^l + (a-bf (A38) 



2i = t + M2 + a - b, q 2 = j + (J.2 - a + b, 



Qz = \- (i2 + ^(i\ + (a + b)\ q a = -- (i 2 - ^A(i\ + (a + bf (A39) 



;i,2 



\ - to) -{a + b)\ 6 = ( J (\ + to) - (a - bf + 2 Ml 



• = j \l[\+to) -(a-6) 2 -2 Ml I . (A40) 



1G 



If W > i.e. Qi > 0, i = 1, 2, 3, 4 then j^l < j and ft > 0, i = 1, 3, hence W is nonseparable for y + (a + 6) 2 > 

3 - A*2 > |a + &| or i < 6 - a - M2- 

Note that the the dimensionality Di given for each item holds for a non-zero choice of the relevant parameters. 
Some eigenvalues Ai may vanish under a special choice of parameters - these subcases are easy to find. There exists 
also symmetric cases 2' and 6' for which the vectors a and b are exchanged. The dimensionality of the local orbits 
remains unchanged, and the formulae for eigenvalues hold, if one exchanges both vectors. 
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